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In models where the breaking of lepton number is spontaneous a massless Goldstone boson, the 
Majoron (J), appears. We calculate the theoretically allowed range for the branching ratios of 
Majoron emitting charged lepton decays, such as Br(/i eJ) and Br(/i eJt), in a supersym- 
metric model with spontaneous breaking of R-parity. Br(/i — > eJ) is maximal in the same region 
of parameter space for which the lightest neutralino decays mainly invisibly. A measurement of 
Br(^ —f eJ) thus potentially provides information on R-parity violation complementary to accelera- 
tor searches. We also briefly discuss existing bounds and prospects for future improvements on the 
Majoron coupling to charged leptons. 



I. INTRODUCTION 

Spontaneous breaking of lepton number leads to a massless Goldstone boson, the Majoron (J) [U, 0, 0]. 
There are two well-known experimental probes for the Majoron: The first is the invisible width of the 
boson, very precisely measured at LEP jj]. The second is neutrinoless double beta decay Q. The NEMO-3 
collaboration, for example, has published limits on half-lives for Majoron-emitting neutrinoless double beta 
decay for a number of isotopes 6]. In addition, there are different astrophysical constraints on the Majoron 
from the cooling of red giant stars and supernovae 0, ^ ■ 

Another interesting possibility to search for Majorons, namely charged lepton decays with Majoron emis- 
sion, has attracted considerably less attention. Indeed, the limits on 1^ IjJ quoted by the Particle Data 
Group H are all based on experimental data which is now more than 20 years old. Probably this apparent 
lack of interest from the experimental side is due to the fact that both, the triplet [2] and the doublet Majoron 
3, are ruled out by LEP data, while the (classical) singlet Majoron model [l| predicts Majoron-neutrino 
and Majoron-charged-lepton couplings which are unmeasurably small. 

Supersymmetry (SUSY) is a theoretically well motivated extension of the standard model. With the first 
data taking of the LHC only months away, searches for SUSY will gain momentum soon. One of the many 
virtues of SUSY is the fact that the minimal supersymmetric extension of the standard model (MSSM) 
provides an interesting candidate for the cold dark matter (CDM), usually assumed to be the lightest 
neutralino (x?) if R-parity is conserved. A stable, electrically neutral lightest supersymmetric particle 
(LSP) will escape the detector and lead to the famous missing momentum signal, upon which standard 
SUSY searches are based. If R-parity is broken, the LSP decays and the CDM candidate is lost. For explicit 
R-parity violation (l^p ), the missing energy signal is degraded, but a larger number of jets and charged 
lepton final states should make discovery of I^p a (comparatively) easy task. In spontaneous I^p {s-l^p ), 
however, the lightest neutralino can decay invisibly through Xi Jv- As pointed out in [l^li if the scale 
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of I^p is low, this decay mode can be easily dominant and s-I^p can be confused with a standard MSSM with 
Rp conserved. 

Here, we revisit U — s- IjJ within spontaneous R-parity violation. Our calculation is based on the model 
of In this model the Majoron is mainly singlet, thus escaping the LEP bounds. This is different from 
the original spontaneous model 3], which used the Icft-sneutrinos to break R-parity. Nevertheless, in our 
model the Majoron can play an important role phcnomenologically. In [l^ U —> IjJ was calculated for a tau 
neutrino mass of rrii,^ ~ MeV. Here we show that (a) despite the fact that current neutrino mass bounds are 
of the order of eV or less, theoretically fj, eJ can be (nearly) arbitrarily large in s-^p , and (b) /i — s- e J is 
large in the same part of SUSY parameter space where the invisible neutralino decay is large, making the 
discovery of R-parity violation at the LHC difficult. Br(/i eJ) thus gives complementary information to 
accelerator experiments. 

At the same time, the MEG experiment has started taking data. MEG is optimised to search for 
Br(/i 67) with a sensitivity of Br(/j, ej) ~ (few) 10^^''. While the impressive statistics of the experiment 
should allow, in principle, to improve the existing bound on Br(/i eJ) 4] by a considerable margin, the 
experimental triggers and cuts make it necessary to resort to a search for the radiative Majoron emission 
mode, Br(/z — > eJj), if one wants to limit (or measure) the Majoron-charged-lepton coupling. We therefore 
also calculate Br(/i — > eJj). 

In the next section, we will briefly discuss s-l^p and give an approximative, analytical estimation of 
Br(Zj IjJ) and Br(Zj — > IjJl)- We then present our numerical results in section (jIII[) . showing the 
correlation of Br(Zj IjJ) with the invisible branching ratio of the lightest neutralino decay. In section 
PVI) we then turn to a brief discussion of existing experimental bounds and comment on how a limit on 
Br(/.t e J7) might be used to put a bound on the Majoron-charged-lepton coupling. We then close with a 
brief conclusion in section (|V|) . 



II. SPONTANEOUS R-PARITY BREAKING 

The model we consider [llj extends the particle spectrum of the MSSM by three additional singlet super- 
fields, 9"^, S and $, with lepton number assignments of L = —1, 1, respectively. The superpotential can be 
written as 

+ KUv^'Hu-h^HdHj' + h^VS+j^''. (1) 

Strictly speaking only v'^ is necessary to spontaneously break Rp. The inclusion of S and $ allows to construct 
a superpotential which purely consists of trilinear terms, thus potentially solving also the ^ problem of the 
MSSM. At low energy various fields acquire vacuum expectation values (vevs). Besides the usual MSSM 
Higgs boson vevs Vd and these are ($) = v^/\/2, {v'^) = Wfl/V2, {S) = vs/V^ and {vi) — wi,./-\/2. Note, 
that Vfi^Q generates effective bilinear terms — hl,vji/^/2 and that vr, vs and ul- violate lepton number 
as well as R-parity. The observed smallness of neutrino masses guarantees that the l^p operators generated 
by these vevs are small. The smallness of implies that either hi, or vji is smalLbut not necessarily both. 

Details of the model, such as mass matrices and couplings, can be found in [ij, [3| , for the phenomenology 
of the LSP decay in this model see 0, [l^ . For brevity in the following we will concentrate on only a few, 
relevant aspects of the phenomenology of this model: The neutrino mass matrix, the lightest neutralino 
decay to Majorons and charged lepton decays. 

The main motivation to study R-parity breaking supersymmetry certainly is that l^p generates neutrino 
masses and thus contains a possible explanation for the observed neutrino oscillation data. For the spon- 
taneous model defined in eq. ([1]), the effective neutrino mass matrix at tree-level can be cast into a very 
simple form 

- i^lZ)i3 = a^^^J + K<^^^j + + ce^ej . (2) 
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Here, = eiVd + VLifJ-, with /i = The coefficients a, h and c are defined as 

m^h?Vch , 1 



AV2Det{MH) 2 



(-hvBVs + -Av^ + hoVdVu), (3) 



DetiMn) 



Det{MH) is the determinant of the (7, 7) matrix of the heavy neutral states (the four MSSM states, B, W 
and Hu,d, plus the three fermionic components of the new singlet superfields of eq. P^) 

Det{MH) = ■^h()h^vl[A{2MiM2^i ~ m^VdVu){-hvRVs + ^Aw^ + h^VdVu) - hom^{vl - v'^f] (4) 

and v"^ — v"^ + v^. The "photino" mass parameter is defined as = g^Mi + g'^M2- Since ^ = 
Det{MH) oc v'^ in the limit of large v^. One can easily fit the observed neutrino masses and angles using eq. 
(H, see [3 

and the short discussion in the next section. 
From a phenomenological point of view the most important difference between spontaneous and explicit 
R-parity violating models is the appearance of the Majoron. The pseudo-scalar sector of the model we 
consider has eight different eigenstates. Two of them are Goldstone bosons. The standard one is eaten by 
the boson, the remaining state is identified with the Majoron. In the limit v^^ ^ vji,vs the Majoron 
profile is given by the simple expression 

(0,0,^,0,^,-^). (5) 



Here, V = + w| and terms of order where v\ — v\., have been neglected. 

Majorons are weakly coupled, thus potentially lead to a decay mode for the lightest neutralino which is 
invisible. Neutralino-Majoron couplings can be calculated from the general coupling X? ~ ~ Pk 

^ = \t{OZ!PL + 0-^^Pu)x',Pl (6) 

Mixing between the neutralinos and the neutrinos then leads to a coupling Xi ~ — J- In the limit 
vrtVs ^ CijfLi one can derive a very simple approximation formula for O^o^^j. It s given by (Iflj 

^ + ^(s'^ii - .9^12) + • • • , (7) 

where the dots stand for higher order terms neglected here and N is the matrix which diagonalizes the 
(MSSM) neutralino mass matrix, e ~ J7j • e and vl — ?7j • vl- Here {Uv)^ is the matrix which diagonalizes 
either the part of the (3, 3) effective neutrino mass matrix, proportional to a or c, depending on which gives 
the larger eigenvalue. Eq. ^ shows that for constant e and wl, Oj^o^^j — > as Wfl goes to infinity. This is as 
expected, since for vu 00 the spontaneous model approaches the explicit bilinear model. We note that, in 
addition to the Majoron there is also a rather light singlet scalar, called the "scalar partner" of the Majoron 
in Sj. The hghtest neutralino has a coupling O^^o^^Sj^ which is of the same order as O^o^^j. Since Sj 
decays to nearly 100 % to two Majorons^his decay mode contributes sizeably to the invisble width of the 
lightest neutralino, for more details see [l(j |. 

The decays k IjJ can be calculated from the general coupling xt ~ Xj ~ Pj!- In the limit of small 
R-parity violating parameters the relevant interaction lagrangian for the U — J — Ij coupling is given by 

^-h{0';^jPL + 0^jPn)hJ (8) 
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with 



Or^j = -^=^[^% + ^(CiA,A,+C'2e.e,+C'3A,e,+C4e.A,)J 

OTI^j = (oZjY- (9) 

The C coefficients are different combinations of MSSM parameters 



„2 

Cl = ^^^i-gWu-VdfJ.^ +VuM2fi) (10) 

^ Det+ 2Det+ 

where Det+ is the determinant of the MSSM chargino mass matrix Z?ei+ ~ ~ ^g^VdVu- Eq. ^ shows 
that one expects large partial widths to Majorons, if vr is low. 

For a charged lepton li, with polarization vector Pi, the decay h — s- IjJ has a differential decay width 
given by 

dcos6' 647rmf ^'•'•^ ^ J ^ * 1' > 

+ \On!,j\\m\ + m| T (m? - m^2)P, cos0) (11) 
+4m,m,i?e(02-^/0^^.,)] 

where is the angle between the polarization vector Pi and the momentum p*,- of the charged lepton in the 
final state, and Pi = \Pi\ is the polarization degree of the decaying charged lepton. 
In the limit nij ~ the expression simplifies to 



dr(Zi I J, J) 



m,. 



since \0'^[j j\'^ oc (ft.^)^ oc m|. The angular distribution of the Majoron emitting lepton decay is thus very 
similar to the standard model muon decay up to corrections of the order {rrij /rrii)^ , which are negligible 
in practice. 

We next consider the decay /i —> eJj. ^ It is induced by the Feyman diagrams shown in fig. ([T]). 




Figure 1: Feynman diagrams for the decay ^ — > eJy. As in the standard model radiative decay /j, ei/vy these 
diagrams contain an infrared divergence for = 0, see text. 



Formulas for the radiative Majoron decays of the r can be found from straightforward replacements. 
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In the approximation ~ the partial decay width for the process /i — s- e J7 can be written as 

r(/i ejj) = ^^\Ol^l^^j\^mf,I{Xmin,ymin) (13) 

where I(xmin ) is a phase space integral given by 

lixrain.Vmin) = J dxdy f {x , y) ^ J dxdy y2^it. X ^-y) ' 
the dimcnsionless parameters a;, y are defined as usual 

a; = , y = ^ (15) 

and Xmin and ymin are the minimal electron and photon energies measured in a given experiment. 

Note that the integral I(xmm,?/mm) diverges for ymin = 0. This infrared divergence is well-known from 
the standard model radiative decay [i evvj, and can be taken care off in the standard way by introducing 
a non-zero photon mass m-y. Note that in the limit me = there also appears a colinear divergence, just as 
in the SM radiative decay. Since in any practical experiment there is a minimum measurable photon energy, 
ymim as well as a minimum measurable photon-electron angle (Oe-y), neither divergence affects us in practice. 
We simply integrate from the minimum value of y up to ymax when estimating the experimental sensitiviy 
of Br(/i — > eJj) on the Majoron coupling. 

In the calculation of the integral T{xmimymin) one has to take into account not only the experimental 
cuts applied to the variables x and y, but also the experimental cut for the angle between the directions of 
electron and photon. This angle is fixed for kinematical reasons to 

cos 6^-1^1 + — . (16) 

xy 



This relation restricts Xmax to be Xmax < 1 as a function of y (and vice versa) and to Xmax < 1 for 
cos > — 1. 

Using the formula for r(/i — )■ eJ), in the approximation We ~ 0, 



r(M-e^) = g:lor;e,7p (17) 



one finds a very simple relation between the two branching ratios 



Br{^ ejj) = —X{Xrain,yrain)Br{pL 6 J). (18) 



We will use eq. p8)) in section pvp when we discuss the relative merits of the two different measurements. 



III. NUMERICAL RESULTS 



All numerical results shown in this section have been obtained using the program package SPheno [If 
extended to include the new singlet superfields P'^, 5* and We always choose the I^p parameters in 



such a way that solar and atmospheric neutrino data [17[ are fitted correctly. The numerical procedure to 
fit neutrino masses is the following. For any random choice of MSSM parameters, we can reproduce the 
"correct" MSSM value of /i for a random value of v^, by appropriate choice of Hq. For any random set of h, 
A, vs and vr, we can then calculate those values of /i^ and v^., using eq. such that the corresponding 
Si and Ai give correct neutrino masses and mixing angles. In the plots shown below we use A.; for the 
atmospheric scale and for the solar scale. 

As shown previously [IQ] if the lightest neutralino is mainly a bino, the decay to Majoron plus neutrino is 
dominant if Vfj is low. This is demonstrated again for a bino LSP in fig. to the left, for a sample point 
using mSugra parameters mo = 280 GeV, mi/2 = 250 GeV, tan/? ~ 10, Aq = —500 GcV and sgn(/i) = +. 
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We stress that this result is independent of the choice of mSugra parameters to a large degree fo!]. A scan 
over V(j, has been performed in this plot, varying in the huge interval [1, 10^] TeV. Large values of lead 
to small values of the constant c in the neutrino mass matrix, see eq. Small c require, for constant 

neutrino masses, large values of e^, which in turn lead to a large invisible width of the neutralino. The largest 
values of (dark areas) therefore lead to the smallest visible neutralino decay branching ratios shown in 

fig. m- 




Figure 2: Branching ratios for visible lightest neutralino decay (left) and branching ratio Br(/i — > ej) (right) versus vr 
in GeV for a number of different choice of v,f, between [1, 10^] TeV indicated by the different colours. Darker colours 
indicate larger in a logarithmic scale. mSugra parameters defined in the text. There is very little dependence on 
the actual mSugra parameters, however, see discussion and fig. ([3]). 

Fig. ([2]), to the right, shows the branching ratio Br(^ eJ) as a function of vji for different values of v^. 
All parameters have been fixed to the same values as shown in the left figure. As the figure demonstrates, 
small values of vr (and large values of w^) lead to large values of Br(/i eJ). This agrees with the analytic 
expectation, compare to eq. 

Our main result is shown in fig. ([3]). In this figure we show Br(/z — > eJ) versus the sum of all branching 
ratios of neutralino decays leading to at least one visible particle in the final state for two different choices 
of mSugra parameters. The similarity of the two plots shows that our result is only weakly dependent on 
the true values of mSugra parameters. We have checked this fact also by repeating the calculation for other 
mSugra points, although we do not show plots here. As expected Br(/i eJ) anticorrelates with the visible 
bino decay branching fraction and thus probes a complementary part in the supersymmetric parameter space. 
An upper bound on Br(/i — > eJ) will constrain the maximum branching ratio for invisible neutralino decay, 
thus probing the part of parameter space where spontaneous R-parity breaking is most easily confused with 
conserved R-parity at accelerators. 

We have checked the points shown in the plots for various phcnomcnological constraints. LEP bounds 
are trivially fulfilled by VLi < vr. Double beta decay bounds on gi„^j ^ are of the order of 10~^ Q and, 
since the coupling g,y^j is suppressed by two powers of R-parity violating parameters, are easily satisfied in 
our model. More interesting is the astrophysical limit on geej- Ref. quotes a bound of geej < 3 • 10~^^. 
Although this bound is derived from the coupling of the Majoron to two electrons, thus constraining actually 
the products v'j^ , and A^, whereas Br(/i eJ) is proportional to Ce^^i and AgA^, it still leads to a (weak) 
constraint on Br(/i eJ), since neutrino physics shows that two leptonic mixing angles are large. This 



^ We will use the symbol g when discussing experimental bounds, to differentiate from the model dependent couplings 
and O"^^ defined in section l lllt . 
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Figure 3: Branching ratios for visible lightest neutralino decay versus branching ratio Br(/i ej) for two mSugra 
points, for various choices of v^, see fig. To the left, same mSugra parameters as fig. ([2J, to the right SPSla'. 

requires that either Ce ^ or Kf. ^ K^. For the case studied in our plots, where ti generate the solar scale, 
tan2 6'0 ~ 1/2 requires eg ~ e^. Numerically we then find that g^ej < 3 • 10 corresponds to an upper 
bound on Br(/K — > eJ) of very roughly Br(/i eJ) < (few) • 10~^. In case the neutrino data is fitted with e 
for the atmospheric scale, the corresponding bound is considerably weaker. 

IV. EXPERIMENTAL CONSTRAINTS AND BR(/i ^ eJ'y) 

The Particle Data Group Q cites [l^ with an upper limit on the branching ratio of Br(/Lt eX°) 
< 2.6 X 10~®, where is a scalar boson called the familon. This constraint does not apply to the Majoron 
we consider here, since it is derived from the decay of polarized muons in a direction opposite to the direction 
of polarization. The authors of [is'l concentrated on this region, since it minimizes events from standard 
model /?-decay. As shown in eq. (jl2p . the Majoron emitting decay has a very similar angular distribution as 
the standard model decay, with the signal approaching zero in the data sample analyzed by [3 ■ Nevertheless, 
from the spin processed data shown in fig. (7) of [1^, which seems to be in good agreement with the SM 
prediction, it should in principle be possible to extract a limit on Br(/.t — > eJ). From this figure we estimate 
very roughly that this limit should be about one order of magnitude less stringent than the one for familon 
decay. For a better estimate a re-analysis of this data, including systematic errors, would be necessary. 

Ref. [l^ searched for Majorons in the decay of tt ^ evJ, deriving a limit of Br(7r evJ) < 4 • 10~^. 
Since the experimental cuts used in this paper [l^ are designed to reduce the standard model background 
from the decay chain tt ^ — > e, the contribution from on-shell muons is reduced by about five orders of 
magnitude. The limit then essentially is a limit on the Majoron-neutrino-neutrino coupling, gu^j, leaving 
only a very weak constraint on the coupling g^ej- Also an analysis searching for Br(/i — > eJ^) has been 
published previously [20]. From a total data sample of 8.15 • 10^^ stopped muons over the live time of the 
experiment [l^l derived a limit on Br(/i c-Jl) of the order of Br(/i eJj) < 1.3 • 10^^. For the cuts 
used in this analysis, we calculate T ~ 10"'^. Thus, see eo. p^ . this limit translates into only a rather weak 
bound Brifj, eJ) < 1.1 • lO'^. 

The MEG experiment [l^ has a muon stopping rate of (0.3 — 1) • 10* per second and expects a total of 
the order of 10^^ muons over the expected live time of the experiment. An analysis of electron only events 
near the endpoint should therefore allow, in principle, to improve the existing limits on Br(/i — > eJ) by an 
estimated (2 — 3) orders of magnitude, if systematic errors can be kept under control. However, the MEG 
experiment, as it is designed to search for Br(/z 67), uses a trigger that requires a photon in the event with 
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a minimum energy of > 45 MeV. MEG data can thus constrain the Majoron-charged-lepton coupHng 

only via searching for T(ii — > eJj). 




Figure 4: The phase space integral for the decay jj, —> eJ^ as a function of a;mi„ for tliree different values of 
J/min = 0.95,0.99,0.995 from top to bottom and for two different values of cosSe^. To the left cos^e^ — —0.99, to 
the right cos 9^^ = -0.99997. 

Fig. ([4]) shows the value of the phase space integral X(a;,„i„, ymin) as a function of Xmin for three different 
values of ymin and for two choices of cos Be-y ■ The MEG proposal describes the cuts used in the search for 
fi as Xjnin ^ 0.995, ymin ^ 0.99 and \tt — 9i,j\ < 8.4 mrad. For these values we find a value of X ~ 

6 • 10"^°. A limit for Br(/i 67) of Br(^ ej) < 10~^^ then translates into a limit of Br(^ eJ) < 0.14, 
obviously not competitive. To improve upon this bound, it is necessary to relax the cuts. For example, 
relaxing the cut on the opening angle to cos^e^ = —0.99, the value of the integral increases by more than 3 
orders of magnitude for Xmin — ymin > 0.95. 

On the other hand, such a change in the analysis is prone to induce background events, which the MEG cuts 
were designed for to avoid. The MEG proposal discusses as the two most important sources of background: 
(a) Prompt events from the standard model radiative decay fj, evPj; and (b) accidental background 
from muon annihilation in flight. For the current experimental setup the accidental background is larger 
than the prompt background. Certainly, a better timing resolution of the experiment would be required 
to reduce this background. For the prompt background we estimate, using the formulas of ^21j], that for a 
total of 10^'^ muon events, one background event from the radiative decay will enter the analysis window for 

^min — ymin 

~ 0.96 for the current cut on cos6'e^. 

A further relaxation of the cuts can lead, in principle, to much larger values for T{xmimymin)- However, 
the search for Br(/Lt e J7) than necessarily is no longer background free. Since all the events from /i eJj 
lie along the line of cosOej defined by eq. (|16p . whereas events from the SM radiative mode fill all of the 
cos space, such a strategy might be advantageous, given a large enough data sample. 

Before closing this section, we mention that tau decays with Majoron emssion are less interesting phe- 
nomenologically for two reasons. First the existing experimental limits are much weaker for taus [2^ 
Br(T fiJ) < 2.3 % and Br(r eJ) < 0.73 %. And, second, although the coupling t — ^ — J is 
larger than the coupling /i — e — J by a factor mr/rrii^i, the total width of the tau is much larger than 
the width of the muon, thus the resulting theoretical predictions for tau branching ratios to Majorons are 
actually smaller than for the muon by a factor of approximately 10^. 
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V. CONCLUSIONS 

We have calculated branching ratios for exotic muon and tau decays involving Majorons in the final state. 
Branching ratios can be measurably large, if the scale of lepton number breaking is low. This result is 
independent of the absolute value of the neutrino mass. The lowest possible values of vp; (at large values of 
v^) are already explored by the existing limit on Br(/i — ^ eJ). 

We have briefly discussed the status of experimental limits. It will not be an easy task to improve the 
current numbers in future experiments. While MEG ^13j certainly has a high number of muon events in the 
detector, a search for Br(/i e J7) instead of Br(/i eJ) suffers from a small value of the available phase 
space integral, given current MEG cuts. An improvement will only be possible, if a dedicated search by the 
experimentalists is carried out. Nevertheless, we believe this is a worthwhile undertaking, since measuring 
a finite value for Br(/i — > e J) will establish that R-parity is broken in a region of SUSY parameter space 
complementary to that probed by accelerator searches. 
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